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Abstract 

Convergence analysis of consensus algorithms is revisited in the light of the Hilbert distance. Tsitsiklis Lyapunov function 
is shown to be the Hilbert distance to consensus in log coordinates. Birkhoff theorem, which proves contraction of the Hilbert 
metric for any positive homogeneous monotone map, provides an early yet general convergence result for consensus algorithms. 
. . , Because Birkhoff theorem holds in arbitrary cones, we extend consensus algorithms to the cone of positive definite matrices. The 

^ I ^ ■ proposed generalization finds applications in the convergence analysis of quantum stochastic maps, which are a generalization 

' of stochastic maps to non-commutative probability spaces. 

■ 
(N : 

L_i ■ I- Introduction 

. The convergence of consensus algorithms in linear spaces has attracted much interest in recent years due to its many 
^5 ' applications in distributed computation and control. Because the matrix defining the update of the algorithm can be given 
. a probabilistic interpretation as a stochastic map or a graph interpretation as an adjacency matrix, the topic has fostered 
ps) ■ a number of interesting connections between control theory, probability theory, and graph theory. Our paper revisits some 
of those connections in the light of a classical fixed point result and uses this link to define a generalization of consensus 
r \ ' algorithms in non-commutative spaces. 

^Z. A key component of our approach is a 1957 theorem of Garrett Birkhoff. Birkhoff theorem provides a general fixed 

point result for homogeneous monotone positive maps defined on closed cones of Banach spaces. The key idea is to use 
f~| . a metric introduced by Hilbert as contraction measure. Hilbert's metric is projective, that is, it measures distances between 

■ homogeneous rays of the cone. 

. In the positive orthant, Hilbert distance between vectors x and y is dH{x,y) = maxlog(a;i/?;i) — mm\og{xi/yi), which 

1 i ' is indeed invariant by scaling of either x or y. Taking y = 1 = [1, . . . , 1]^, this provides a natural distance to consensus; 

in probabilistic terms this corresponds to a uniform density vector In fact, Hilbert distance to consensus coincides (in log 
coordinates) with the Lyapunov function introduced by Tsitsiklis to study the convergence of consensus algorithms. Row- 
stochastic matrices define linear positive maps in the positive orthant. As a consequence, Birkhoff theorem provides an early 
convergence result for consensus algorithms. 
\^ , Because Birkhoff theorem applies in arbitrary cones, it opens the way to several generalizations of consensus theory. The 

■ present paper focuses on the cone of hermitian positive definite matrices. Our motivation stems from the generalization of 
. probability theory in the quantum (non-commutative) setting. While stochastic maps operate on probability vectors (elements 

■ of the positive orthant), quantum stochastic maps operate on probability matrices (elements of the cone of hermitian positive 
, semidefinite matrices). Quantum stochastic maps are unital maps in this cone, i.e. they map identity on the identity. Their 

\ I ' dual (Kraus maps) are trace-preserving and appear frequently in the study of open quantum systems. 
I> . In the cone of hermitian positive definite matrices, Hilbert distance from a matrix X to the identity matrix is dn {X, I) = 
' log(Amax) — log(Amin)- As in the positive orthant, this distance is projective, i.e. invariant by scaling (of either X or /). For 
unital maps, it provides a natural distance to the uniform probability density matrix /, which defines "consensus" in this 
context. Because quantum stochastic maps define positive linear monotone maps on the cone of positive definite matrices, 
Birkhoff theorem shows that the distance to consensus can only decrease when measured with the Hilbert metric. The paper 
explores some implications of this basic property in the context of quantum channels. The preliminary examples discussed 
in this paper indicate that the analogy between the classical consensus theory in commutative spaces and the proposed 
consensus theory in non-commutative spaces offers new avenues to characterize the convergence properties of quantum 
channels. In particular, Tsitsiklis approach to consensus theory proved especially useful in asymmetric and/or time-varying 
settings where a more conventional approach based on quadratic Lyapunov functions fails. Likewise, the Hilbert distance as 
a measure of distance to consensus in a quantum setting could prove particularly useful in the context of asymmetric and/or 
time-varying quantum maps. 

The paper is organized as follows: the convergence theory of consensus algorithms is briefly revisited in Section 2. Birkhoff 
theorem is recalled in Section 3. Its interpretation in the framework of consensus algorithms is discussed in Section 4. Section 
5 extends the framework to the cone of positive definite matrices. The notion of non-commutative consensus, motivated by 
quantum applications, is introduced in Section 6. Section 7 contains concluding remarks and future perspectives. 
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II. Consensus algorithms and Tsitsiklis Lyapunov function 
Linear consensus algorithms give rise to time-varying systems 

x{t + l) = A{t)x{t), x(t)€R" (1) 

where for each t = 0,1, . . . the matrix A{t) is row-stochastic, i.e. its elements are non-negative (aij(t) > 0), and it preserves 
the unit element 1 = (1, . . . , 1)^ (A{t)l = 1). The graph interpretation is that N nodes exchange information about a scalar 
quantity Xi{t) along communication edges weighted by the coefficients aij{t). {aij{t) = means that no information is 
transferred from node j to node i at time t). The update of a given node is a weighted average of its own value and the 
values communicated by neighbors: 

n n 

x,{t + l) = ^a,j{t)xj{t), ^a,j{t) = l (2) 

Remarkably, uniform convergence of the time-varying system ([T]i can be characterized quite precisely. The early analysis 
of Tsitsiklis [9] rests on the basic but fundamental observation that the (time-invariant) Lyapunov function 

V{x) = max Xi — min Xi (3) 

l<?'<n l<?'<n 

is non-increasing along the solutions of ([T). The proof that it decreases strictly over a uniform horizon under appropriate 
assumptions only involves elementary calculations (see e.g. [4] for details). 

An example of general convergence result for ([T]i is as follows: 
Theorem 1: Let A{t), t > be a sequence of n x n row-stochastic matrices and suppose that the following conditions 
hold: 

• all nonzero Uij (t) are larger than a certain a > 0. 

• All diagonal elements satisfy an (t) > 

Then the solution of Q converges exponentially to an equilibrium co-linear to 1 if and only if there exists a finite horizon 
T > such that for all to e N, there exists a node connected to all other nodes across [to, to + T]. 

The proof of (variants) of this result has appeared in various papers, including important contributions by Tsitsiklis [9], 
Jadbabaie et al. [1], and Moreau [10]. 

It is also remarkable that the convergence analysis of the linear system ^ cannot be established by means of a time- 
invariant quadratic Lyapunov function. The recent paper [13] provides the explicit construction of a sequence of eight 
matrices A{t) that satisfies the assumptions of Theorem 1 but that does not admit a common time-invariant quadratic 
Lyapunov function. In other words, the non-quadratic nature of the Lyapunov function ^ is essential and a distinctive 
feature of the analysis of consensus algorithms. 

III. BiRKHOFF Theorem 

As will be shown in the next section, the Lyapunov function (O is closely related to Hilbert metric. This connection 
follows from a result of Birkhoff [3] that we now summarize, following the exposition in [5]. 

Let A" be a real Banach space and let /C be a closed solid cone in X, that is, a closed subset K, with the properties: (i) 
intAC, the interior of /C, is not empty, (ii) /C + /C C /C, (iii) A/C C /C for all A > 0, and (iv) K, D — /C = {0}. The partial 
order x <y means that y ~ x ^ JC. 

For x,y G ICq = /C\{0}, define M{x,y) = inf{A x — \y < 0} and m(x,y) = sup{A : x — Xy >z 0}- Then Hilbert 
metric d{-,-) is defined in ICq by 

d{x, y) = \og{M{x, y)/m{x, y)} (4) 

Because of the invariance property d{ax,Py) = d{x,y) for all a,/3 > 0, Hilbert metric is, strictly speaking, a projective 
metric, i.e. a distance between the equivalence classes (rays) [x] and [y], where [x] = {ax | a > 0}. 

Birkhoff theorem is a contraction result for positive maps in JC. A map A : JC ^ K- m said to be non-negative; a map 
A : int/C int/C is said to be positive. If A is positive, the projective diameter A(A) of A is defined by 

A(yl) = s.\vp{d{Ax, Ay) : x,y e int/C} . 

Theorem 2 [Birkhoff, 1957]: Let A a map in K. that satisfies the following properties: 

(i) A is positive, i.e. A : int/C — > int/C. 

(ii) A is homogeneous of degree p in int/C, i.e. A{Xx) = X^A{x) for all A > 0. 

(iii) A is monotone, i.e. x ^ y => A{x) ^ A{y). 

Then A is a contraction for the Hilbert metric and the following holds: Vx,y G int/C : 
d(A(x),A(y)) <tanh(iA(A))d(.T,y) 



It should be noted that Banach contraction map theorem applied to Birkhoff theorem provides a far-reaching generalization 
of Perron-Frobenius theorem: if A(^) < oo and if the metric projective space E = (int/C\"', d) is complete, then there exists 
a unique eigenvector of A in E. 

Hilbert metric is only one among the many metrics contracted by positive monotone maps. Another such (closely related) 
metric is Thompson metric dT{x,y) ~ \ogma.x{M{x,y),m~^{x,y)}. However, Hilbert metric typically provides the best 
contraction ratio, see [5] for details. 



IV. Birkhoff Theorem in the positive orthant and Consensus algorithms 

The most direct application of Birkhoff theorem is in the positive orthant. In this case, X = R" and K, = {(2:1, ■ ■ ■ ,Xn) ■ 
Xi > 0, I <i < n}. Then M{x,y) = maxi{xi/yi) and m{x,y) = mini{xi/yi) and hence Hilbert metric expresses as 

d{x,y)^ log . 7 ■ 

By definition, stochastic matrices A{t) define positive monotone maps in /C. By linearity, they are also homogeneous of 
degree one. 

Because 1 is a fixed point for all A{t), Birkhoff theorem implies that d{x, 1) < tanh(iA(A(t))) d{A{t)x, 1) for all t. 
This means that for each x e int/C, Birkhoff theorem provides the Lyapunov function 

Vb{x) = d{x, 1) = log ^ maxlog(a;i) - minlog(2;i) (5) 

This is exactly the Lyapunov function ^ in log coordinates. Both Lyapunov functions provide a measure of the projective 
distance from [x] to [1]. Tsitsiklis Lyapunov function is translation-invariant {V{x + Al) = V{x)) while Birkhoff Lyapunov 
function is scaling-invariant {Vb{Xx) = V{x) for A > 0). They both provide a measure of the diameter of the convex hull 
of (.Ti, . . . , Xn), which is the Lyapunov function proposed by Moreau [10]. 

Birkhoff theorem is thus an early contribution that points to two important facts: (i) the natural distance to study the 
convergence of consensus algorihms is not quadratic and (ii) the convergence analysis should concentrate on the contraction 
of the interval [xmin, a;max]- 

Birkhoff theorem also provides a contraction coefficient. It is given by the diameter of A{t) 

A(A(t)) = sup(log(^^47H7T) :1<*'J'P'9<" 

When the diameter is finite, the contraction coefficient is strictly smaller than one. Uniform convergence only requires 
contraction over a uniform time-horizon. As a consequence, if A{t) does not have a finite diameter at each instant, it is 
sufficient for uniform convergence that the finite product A{t) = A{t + T) . . . A{t + l)A{t) has a finite diameter 

In the symmetric case A{t) = A'^{t), the existence of a finite horizon over which the diameter is finite is closely linked 
to the existence of a finite horizon over which the communication graph is connected. This link is studied in detail in [1]. 
It provides a close connection between Birkhoff theorem and Theorem 1 in the symmetric setting. Symmetric stochastic 
matrices are called doubly-stochastic matrices (they are row-stochastic and column-stochastic) and symmetric graphs appear 
most notably when the edges are undirected. 

In contrast, the asymmetric setting (that occurs in directed communication graphs) includes situations where, despite a 
uniform convergence to consensus, the diameter is infinite over any finite horizon. In those situations, illustrated in the next 
example, Birkhoff theorem does not provide a strict contraction measure even though elementary explicit calculations show 
the contraction of [xmin, a^max] under the assumptions of Theorem 1. 
Example 1: A consensus map with infinite diameter. The matrix 

A= ^, ) with7e(0,l). 

satisfies the assumptions of Theorem 1. For an initial condition (a;i, X2) in the positive orthant, the asymptotic (consensus) 
value is {xi,xi). Nevertheless, the diameter of A is infinite. Because A is lower-triangular, any finite power A'^ also has 
an infinite diameter □ 
Quantifying the convergence rate with other coefficients than Birkhoff diameter remains a topic of research, see e.g. [2] 
and references therein. 

It should be emphasized that Birkhoff theorem is not restriced to linear maps (see [7], [10] for examples of nonlinear 
consensus algorithms). 



V. BiRKHOFF Theorem in the cone of positive definite matrices 

We now turn to the application of Birkhoff theorem in the cone of positive definite matrices. In this case, X = {X = 
X* £ is the set of hermitian matrices, where * denotes complex conjugate transpose, and K, = {X ^ | X e X}. 

Then for X,Y yO, AI{X,Y) = Xmax{XY^^) and m{x,y) ~ Xi^in{XY^^). Hence Hilbert metric expresses as 

The reader will note that an equivalent expression is 

which shows the close connection of Hilbert metric to the natural invariant metric 

dRiX,Y) = \\logiY-'/^XY-'/')\\F^ZLi^ogK{Y-'^'XY-'^'). (7) 

In contrast to the Hilbert metric, the metric (|7| is a Riemannian metric. Both metrics share congruence invariance properties: 
d{FX F* , FYF*) = d{X,Y) for any invertible linear map F in C"^". The Riemannian metric coincides with the Fisher- 
information metric in statistics [14] and with the self-concordant logarithmically homogeneous barrier (— logdetA) in 
optimization [11]. It makes the cone of positive definite matrices geodesically complete. 

Birkhoff theorem implies that any homogeneous monotone positive map A on the cone of positive definite matrices 
contracts the metric (|6]). 



VI. Non-commutative consensus 

A. Generalizing stochastic maps to non-commutative spaces 

A key feature of consensus algorithm ([T]) on K" is the stochasticity of maps A{t), i.e. (a) they have non-negative elements 
and (b) A{t)\ — 1. An interesting analog of stochastic maps on the cone JC of positive semi-definite hermitian matrices of 
order n is given by the dual maps of Kraus maps describing quantum channels (see [12] chapter 8 for a tutorial presentation 
of quantum channels; see [8] chapter 4 for a physical presentation of decoherence and Kraus maps). These dual maps are 
unital completely positive maps, which admit the following characterization [6]: a linear map $ on /C is unital completely 
positive if and only if it admits the expression 

VXe/C, = with ^V*V,=I (8) 

i i 

for matrices Vi E C"^". Completely positive maps, characterized by the first line of (O, map K, into itself. They are 
positive linear maps on K, in the terminology of the present paper Indeed, X y implies that, for w € C", w*^{X)w = 
Y,, {V^w)*X{Viw) = if and only if V^w = Vi But then {V^w)*V^w = and thus w = since ^^^^ = ^■ 
Complete positivity of $ generalizes the non-negative elements property of stochastic matrices. The additional unital property 
$(/) = / is the analog of A{t)\ = 1 for stochastic matrices. 

This suggests the following non-commutative extension of consensus algorithms: 

X(t + 1) = X(0) = Xo e int/C (9) 

where for each t, the map $t is of the form dHJ. The linear consensus algorithm on M" can be obtained as a special case 
of dill, see Section lyTcl 

Note that, in the same way as algorithm ([T]i can be extended from the positive orthant to the entire space M", the non- 
commutative generalization (|9]l can be extended to arbitary hermitian matrices in X ~ {X = X* e C"^"}. This is a 
consequence of the translation invariance property of the algorithm: the property $t(/) = / and linearity imply that, if X{t) 
is the solution of ^ with X{0) — Xq e /C, then the solution of (|9]l with initial condition Xq + al equals X{t) + a I, for 
any a G M (in particular a < 0). 



B. Birkhoff theorem and convergence of non-commutative consensus 

By definition, $t is a linear monotone positive map on the cone of positive definite matrices (i.e. satisfying conditions (i)- 
(iii) of Theorem 2). The following generalizes to the non-commutative setting (|9]l the basic convexity property, {^^(t+l), i = 
l,2...Ti} C [max{xi{t)), mm{xi{t))] for solutions of ([T]i. 

Lemma 1: Let ^ be a map of the form (HJ. Then Ai„i,i($(X)) > X^-ainiX) and Ainax('I'(^)) < ^maxiX) for any X = 
X* e C"^". 

Proof: Define Yi := X - {Xrain{X) - e)I, for any < £ < 1. Then Yi e int/C implies $(Yi) G int/C. Since 
^{Yi) := — (Aniin(X) — e) / we have the first bound when e — > 0. The second bound is obtained analogously with 

Y2:=-X + iX^UX)+e)L □ 



The following result is a direct consequence of Birkhoff theorem and applying Lemma 1 recursively at each t. 
Theorem 3: Let <l>t, t >0 be a sequence of maps of the form (jS]). Then the Lyapunov function 

is non-increasing along the solutions of Denoting $[( 4+7^] := '^t+T o • • • o $4+1 o $j, the Lyapunov function satisfies 

V{X{t + T)) < tanh(iA($[t,t+T])) V'(^) 

In particular, if the maps ^[t^j^x] have a finite projective diameter for some finite horizon T, then the decrease ofV is strict 
over T iterations of and the solution of ^ uniformly converges to a point in the set {A/ | A £ [Aniin(Xo): Amax(Xo)]}- 

The following proposition is useful to characterize the diameter of the application. 

Proposition 1: Let $ an abbreviated notation for the map $[( 4+7-] corresponding to a finite number of iterates of ^ as in 
Theorem 3. 

(a) Define i?(l>) = sup{log ^"""[l^^jj I X e int/C}. Then i?($) < A($) < 2 

(b) If ^(P) has full rank for all rank-1 projectors P V = {xx* \ x € C", ||a;||2 = 1}, then the superior value in i?($) 
is obtained for X Cz V and R{^) ~ R < 00. 

Proof: (a) A($) > i?($) follows by choosing F = / in the definition of A. A(l>) < 2i?($) follows from the triangle 
inequality on the Hilbert metric. 

(b) If $(P) has full rank for all P e V, then i?p($) := sup{log ^"'°"(||p|| I P e iiitT'} takes a finite value, since V 

is compact and P 1-^ log ^"""If |pj| is continuous. Every X <E K. can be expressed as a sum of positively weighted rank-1 
projectors. Since Ainax(-''^) and Aniin(-''^) are respectively convex and concave functions of X e /C and d{X,I) is invariant 
by scaHng of X, we directly have i?p(l>) = R{^)- □ 

Note that a finite diameter for the finitely iterated map is only a sufficient condition to prove strict contraction. Very much 
as in the classical consensus setting, there are situations where the strict decay of the Lyapunov function will be established 
by direct computations despite an infinite diameter over any finite horizon. 

C. Consensus in as a special case of non-commutative consensus 

There is an obvious bijection between the positive orthant of M" and the set of diagonal positive semidefinite matrices 
V = {X £ R"^" \ X >z 0, Xij = if i 7^ j}. In this section we show how consensus in M" can be recovered as a 
special case of the non-commutative generalization, when the state is restricted to diagonal matrices in some fixed basis. 
These diagonal matrices do all commute with each other; this justifies the appellation "non-commutative" consensus for the 
generalized case of non-diagonal matrices. 

Denote by S the set of all permutation matrices, i.e. matrices with exactly one element equal to 1 on each row and 
column, and all other entries 0. Note that S*S = I for any S £ S. Hence a map 

X <i>s{X)^Y^ Wi S*XS^W^ with S^eS, W^eVyi md^W^ =1 (10) 

i i 

is of the form dHJ. Moreover, $5 defined by ( fTOl i maps T> into itself. Indeed, the map X S*XSi permutes in the same 
way both the rows and columns of X; for X £ T> this amounts to permuting its diagonal elements. The Wi further scale 
each element. This suggests a way to write the classical linear consensus dU in the form of (|9]l on 2? where Xi of ([T]i is 
identified with xa of (|9]l. For instance, take Si to be the permutation sending Xjj to position (j © i, j © i) for all j, where 
ffi denotes addition modulo n. Further take element (j © i,j ffi i) of Wi equal to y/aj(si, i- Then ( fTOb with this setting is 
equal to ([T]); in particular, i^Djj ~ J2i ^ji — 1- (This construction is not unique.) 



D. Birkhoff theorem on dual and convergence of iterated Kraus maps 

For simplicity, we consider a fixed (time-invariant) map $ for the rest of the paper The adaptation to the time-varying 
case is straightforward but left for an extended version of the paper 

The dual map 'I' of $ defined by ([8]) is the usual framework to describe the evolution of the density matrix characterizing 
a state in a quantum channel, see e.g. [12]. It is called a Kraus map and it admits the form 

vze/c, -^{z) = ^ViZV* (11) 



with still ^ j V*Vi = /. The unital property of $ becomes a trace-preserving property for the dual map Kraus maps are 
thus trace preserving completely positive maps. In quantum applications, the density matrix typically evolves according to 

Z{t + l) = '^{Z{t)), Z{t) e ICi := {X etc \ti{X)^l} (12) 

where the Vi defining are matrices in C"^". From Theorem 3 and Proposition 1, we can deduce the following convergence 
result for 4'. 

Theorem 4: Let 5' defined by Allh Assume that N -iteration of its dual map $ defined by ^ admits a finite projective 
diameter A($^) < oo, for some integer N . Then, there exists a unique fixed point Z of ^ in ICi, and for any Zq G JCi, we 
have limi_>+oo Zit) = Z. For any initial Xt) G X of Xit) defined by (|9|l, we have limt_>.+oo ^(^) = tr(ZXo)/. Moreover 
the convergence is exponential. 

Proof: $t satisfies assumptions of Theorem 3. Thus for any initial hermitian matrix Xq, X{t) converges exponentially 
towards f{Xo)I where f{Xo) G [X^in{Xo), Aniax(-'^o)] is the consensus value, which depends on Xq. Since <bt is a linear 
map on X, f is also a linear map on X. Therefore there exists a unique hermitian matrix Z such that f{Xo) = tr(ZXo) 
WXq e X. For all t we have, by duaUty, tr(Z(t)Xo) = tr(ZoX(t)). Thus for any Xq € X, tr(Z(t)Xo) converges 
exponentially towards tr(Zo f{Xo)I) = f{Xo) tr(Zo) = f(Xo) = tr{ZXo). Thus Z{t) converges exponentially towards Z, 
which is positive semidefinite since K, is closed. ti{Z) = 1 since ti{Z{t)) = 1. □ 

The relation between primal and dual maps can also be made for the classical consensus in M". Indeed, viewing (T^ as 
the dual, the primal is 

z(t+ 1) = A* zit), zit) e M" (13) 

where A* is now column stochastic. This implies that J^i is conserved. Analogously to Theorem 4, the scalar z*Ax 
can be expressed as evolving either on the primal or on the dual; thus if ([T]i converges to consensus for all ccg, then (fT3] l 
converges to some z (not necessarily consensus). 



E. Illustration on quantum channel applications 

Example 2: spin rotations. We consider an open quantum channel where a stochastic environment influences a system 
with density matrix p e C^^^ — that is a two-level system like e.g. a spin or qbit (see [8], [12]). Generalizing [12], we 
assume that at each iteration, two different "spin-rotations" can be applied to the system with probabilities p and 1 — p 
depending on the stochastic environment. The expectation Z E C^^^ of all stochastic (Monte-Carlo) trajectories for p is 
then characterized by 

Zit + 1) = VoZVo* + ViZV{ (14) 

■t- e"^ \ f cos(/3) i sin(^) \ 

with Vo:=r. = Vp[ e"'" j ' = [ ^ sin(% cos(/3) ) ' ^''^ 

Note that J^k ^k^k = X^A: Vk = ^- This implies that the quantum stochastic map is both trace-preserving and unital. 
In fact the Kraus map $ and its dual ^I* are equal, modulo inversion of the angles a and /3. 

Therefore, we can view ( fT4] i directly (instead of its dual) as a non-commutative consensus iteration $. We compute that 
A($'^') < oo for fc > 2, except if (a) a ~ OmodTr or (b) /3 ~ OmodTr, or (c) (2a, 2/3) = (0, 0)mod7r. Thus except for these 
special cases, we can apply Proposition 1 and Theorem 3 to conclude that Z converges to ^ /. The physical interpretation 
is that, when stochastically alternating between two generic rotations, the probability distribution for the state becomes 
completely uniform after many iterations. 

A Kraus map that is both unital and trace-preserving can be considered as the analog of a doubly stochastic map in 
classical consensus theory. The consensus value 1/2 for a trace-perserving two-dimensional quantum map is reminiscent of 
the average-preserving property of doubly stochastic maps. 

For the particular value /3 = f , the set of diagonal matrices X is invariant under ( fT4l i. In fact. 



pxii + (1 - p)a;22 

PX22 + il - p)xii 



for a diagonal matrix X = diagj^n, X22} which can be interpreted as a classical consensus algorithm for the doubly- 
stochastic map 



A 



P (1-p) 
il-p) p 



□ 



Example 3: spontaneous emission. Consider a two-level quantum system as in Example 2. The expectation trajectory of 
the system undergoing the stochastic process of spontaneous emission is described by 

V - V -( ^ J ^ V - V _ 7 

*Q • — no emission y \/^ / ' ' Photon emission 

with small 7 equal to the quotient of the chosen time discretization step by the expected lifetime of the unstable excited 
state. The dual of the Kraus map is 



xii X12 \ ^fluv\-l ^11 Xi2\/1- 7^ 



X22 ) V Xl2\jl-l'^ X22{l-1^)+Xiir 

From the simple form of the map, it is easy to guess asymptotic convergence to the consensus state xnl. This property is 
reminiscent of the consensus property of a directed graph where the root agent eventually imposes its value to all others. 
Indeed, restricting <i> to diagonal X results in the commutative consensus map of Example 1. (Following the construction 
proposed in Section VI.C, one has Vq ~ SqWq and Vi — SiWi with Wo = diag{0,7} and Wi = diag{l, \/l ~ 7^}-) 

As in Example 1, this is a situation where the diameter of the quantum stochastic map is infinite for any finite iteration, 
yet the contraction of spectrum range can be computed explicitly. The calculation gives 

[A„,in(3>(X)), A„,ax($(^))] - [XminiX) + p+, A„,ax(^) - P-] 

for p± = 7^ (^£11.^^22 -J 2 _|_ |2;i2p ± ( ^"~^^^ )] . This shows the strict decay of the Lyapunov function V at each iteration. 

□ 



VII. Conclusion 

We studied some implications of Birkhoff theorem for the convergence analysis of consensus algorithms in two different 
cones; the positive orthant of the eucUdean space and the cone of hermitian positive definite matrices. In the positive orthant, 
the application of Birkhoff theorem to linear stochastic maps provides an interpretation of Tsitsiklis Lyapunov function as 
the Hilbert distance to consensus in log coordinates. This connection does not seem to have been made in the literature. In 
the cone of positive definite matrices, the application of Birkhoff theorem provides a framework to study the convergence of 
quantum stochastic maps. This framework suggests to study the contraction of the interval [Amin, A,nax] over finite horizons. 
To the authors' knowledge, this approach is novel. The analogy with the classical consensus problem suggests that the 
potential of this approach is in situations where the stochastic map is not self-adjoint (the analog of directed communication 
graph in classical consensus theory) and/or in situations where the stochastic map is time-dependent (the analog of time- 
varying communication graphs in the classical consensus theory). These are the situations where common quadratic Lyapunov 
functions may fail to exist. The analog of Theorem 1 in the non-commutative setting remains to be characterized. 
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